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We look at gravitational attraction in simplicial gravity using the dynamical triangulation method. On the 
dynamical triangulation configurations we measure quenched propagators of a free massive scalar field. The 
masses measured from these propagators show that gravitational attraction is present. 



1. Introduction 

One of the approaches to quantum gravity that 
has recently raised increasing interest is that of 
dynamical triangulation. More details on this 
method by itself can be found in e.g. [1,2]. 

The one thing that everybody knows about 
gravity is that it is an attractive force. This led 
us to the wish of actually seeing some of this at- 
traction in dynamical triangulation. 

2. Description of the model 

We look at the behaviour of a free scalar field 
(f) with bare mass mo in a quantum gravity back- 
ground. The euclidean action of this system in 
continuum language would be 
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where Aq is the bare cosmological constant, R is 
the scalar curvature and Go is the bare Newton 
constant. 

We take (f) as a test particle here, i.e. the back 
reaction of the field (f) on the metric is not taken 
into account. It is seen in other simulations [3] 
that such matter has little infiuence on the gravity 
sector of the theory. 



We will use the following notation for expecta- 
tion values of an observable A. 



IV^Aexpi-S[g,^]) 
fV<j,eM-S[g,<i>]) ' 
fVgAeM-S[g]) 



JVg exp{-S[g]) ' 
The quenched expectation value is then 
{A) = {{A)^),. 
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We can now look at propagators in a fixed ge- 
ometry. The one particle propagator, denoted by 
G(x), is defined as 



G(x) = {c^^c^ 



(7) 



where is an arbitrary point. The connected two 
particle propagator will then be the square of the 
one particle propagator 
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Letting the metric fiuctuate, we take the av- 
erage of the propagators over the different met- 
rics. Because of reparametrization invariance, the 
average (G{x)) g will not depend on the place x. 
Therefore, we look at averages at fixed geodesic 
distance d 



{G{d))g = \ IVg eM-S[g]) x 



d°xG{x)b{d{x,Q) - d), 



(9) 
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where d(x, y) is the minimal geodesic distance be- 
tween X and y. For a massive particle, we expect 
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Figure 1. The two particle and the square of the 
one particle propagator versus the geodesic dis- 
tance for three different bare masses in two di- 
mensions. The vertical scale is logarithmic. 

this propagator to fall off exponentially as 

{G{d))g (X rf«exp(-mrf), (10) 
{G{d))] (X rf2«exp(-2mrf), (11) 

with some power a and the measured mass m. 
The two particle propagator will behave similarly 
as 

{G{df)g oc d^eM-Ecd), (12) 

where is the energy of the two particle com- 
pound. It this energy turns out to be less than 
two times the mass of a single particle, the dif- 
ference can be interpreted as a binding energy 
between the particles. This would show gravita- 
tional attraction between them. 

3. Implementation 

We have run numerical simulations with both 
two and four dimensional dynamical triangula- 
tions. 

In two dimensions the volume can be kept con- 
stant, and for fixed topology no parameters will 
be left. We used systems of 32k and 64k triangles 
with the topology of the two-sphere. 



Figure 2. The effective binding energy E}, as a 
function of the geodesic distance for three differ- 
ent bare masses in two dimensions. 

In four dimensions the analogue of the contin- 
uum gravitational action is 

^[^'] = T77^ I d^^Vai^^o-R), (13) 

iDTTCro J 

K4N4-K2N2, (14) 

where N2 and N4 are the number of triangles and 
four-simplices repsectively. We used systems of 
about 16k simplices and the topology of the four- 
sphere. To keep the number of simplices around 
the desired value, we varied K4, increasing it when 
the volume became too large and vice versa. 

The parameter K2 is inversely proportional to 
the bare gravitational constant Go- As in other 
work [1-4], we see indications of a second order 
phase transition as K2 varies. 

On each dynamical triangulation configuration 
we then calculated the propagator 

G(x) = (D^ + m^)-; (15) 

of the scalar field, using the algebraic multigrid 
routine AMG1R5. 

Eventually we hope to be able to extract the 
renormalized Newton's constant Gr at the criti- 
cal value of K2, e.g. according to the nonrelativis- 
tic formula 

Ei = ^Glm'. (16) 
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Figure 3. As in figure 1, but in four dimensions. 
K2 = 1.2, which corresponds to the weak coupling 
(elongated) phase. 



Figure 4. As in figure 2, but in four dimensions. 
Again, K2 = 1.2. 



(12) we have 



4. Results 

In figure 1 we see the results in two dimensions 
for three different bare masses. Each pair of lines 
corresponds to one bare mass. In each pair the 
upper line is ln{G(d)^) (the two particle propaga- 
tor) and the lower line is ln{G(d))^ (the product 
of two single particle propagators). 

Their is clearly a difference in slope between 
the lines in each pair. This shows that the en- 
ergy of the two particle compound is less than 
two times the mass of a single particle and con- 
sequently that there is a positive binding energy 
between the particles. 

Figure 3 shows similar data in four dimensions. 
The propagators at the lowest mass show large 
fiuctuations, but at the higher masses there is 
again a clear difference in slope. 

In the two dimensional case, one might expect 
not to see any attraction, because of the absence 
of dynamics in the classical system. In the quan- 
tum case, however, a non-trivial theory can be 
seen due to the conformal anomaly [5] . 

From these data, we can now calculate the 
binding energy of the particles. From (11) and 
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Figure 2 shows this quantity as a function of 
the geodesic distance. The three curves again 
correspond to the three different bare masses. Al- 
though the result is not yet very accurate, it is 
clear (and fits to these curves show) that the bind- 
ing energy goes to a non-zero value. 

Figure 4 shows the same data in the four di- 
mensional simulations. Obviously, these data are 
very preliminary and we need more statistics. 
Nevertheless, this figure also indicates a non-zero 
binding energy. 
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